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We present a model-independent method to test for scale-dependent non-Gaussianities in combi- 
nation with scaling indices as test statistics. Therefore, surrogate data sets are generated, in which 
the power spectrum of the original data is preserved, while the higher order correlations are partly 
randomised by applying a scale-dependent shuffling procedure to the Fourier phases. We apply this 
method to the WMAP data of the cosmic microwave background (CMB) and find signatures for 
non-Gaussianities on large scales. Further tests are required to elucidate the origin of the detected 
anomalies. 
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Inflationary models of the very early universe have 
proved to be in very good agreement with the observa- 
tions of the linear correlations of the cosmic microwave 
background (CMB). While the simplest, single field, 
slow-roll inflation [U [21 [3] predicts that the temperature 
fluctuations of the CMB correspond to a (nearly) Gaus- 
sian, homogeneous and isotropic random field, more com- 
plex models may give rise to non-Gaussianity [H |5l El [7] . 
Models in which the Lagrangian is a general function of 
the inflaton and powers of its first derivative [HI [9^ can 
lead to scale-dependent non-Gaussianities, if the sound 
speed varies during inflation. Similarily, string theory 
models that give rise to large non-Gaussianity have a 
natural scale dependence [10 . If the scale dependence 
of non- Gaussian signatures plays an important role in 
theory, the conventional (global) parametrisation of non- 
Gaussianity via fjsfL is no longer sufficient to describe 
the level of non-Gaussianity and to discriminate between 
different models, /nl must at least become scale depen- 
dent - if this parametrisation is sufficient at all. But first 
of all such scale-dependent signatures have to be identi- 
fied. 

Possible deviations from Gaussianity have been investi- 
gated in studies based on e.g. the WMAP data of the 
CMB (see [TT] and references therein) and claims for 
the detection of non-Gaussianities and other anomalies 
(see e.g. [Il[T3[l4l[T5l[Ii[ni[Tl[ll[20]) have been 
made. These studies have in common that the level of 
non-Gaussianity is assessed by comparing the results for 
the measured data with a set of simulated CMB-maps 
which were generated on the basis of the standard cos- 
mological model and/or specific assumptions about the 
nature of the non-Gaussianities. 

On the other hand, it is possible to develop model- 
independent tests for higher order correlations (HOCs) 
by applying the ideas of constrained randomisation [21, 
[22l [23] , which have been developed in the field of non- 
linear time series analysis [24|. The basic formalism is 
to compute statistics sensitive to HOCs for the original 



data set and for an ensemble of surrogate data sets, which 
mimic the linear properties of the original data. If the 
computed measure for the original data is significantly 
different from the values obtained for the set of surro- 
gates, one can infer that the data contain HOCs. 
Based on these ideas we present in this Letter a new 
method for generating surrogates allowing for probing 
scale-dependent non-Gaussianities. Our study is based 
on the WMAP data of the CMB. Since our method in its 
present form requires full sky coverage to ensure the or- 
thogonality of the set of basis functions Yim we used the 
five-year " foreground- cleaned" Internal Linear Combina- 
tion (ILC) map (WMAP5) [2F generated and provided^ 
by the WMAP-team. For comparison we also included 
the maps produced by Tegmark et al. [26|[22 , namely the 
three year cleaned map (T0Hc3) and the Wiener-filtered 
cleaned map (T0Hw3)^, which were generated pursuing 
a different approach for foreground cleaning. Since the 
Gaussianity of the temperature distribution and the ran- 
domness of the set of Fourier phases are a necessary pre- 
requisite for the application of our method we performed 
the following preprocessing steps. First, the maps were 
remapped onto a Gaussian distribution in a rank-ordered 
way. By applying this remapping we automatically focus 
on HOCs induced by the spatial correlations in the data 
while excluding any effects coming from deviations of the 
temperature distribution from a Gaussian one. 
To ensure the randomness of the set of Fourier phases we 
performed a rank-ordered remapping of the phases onto 
a set of uniformly distributed ones followed by an inverse 
Fourier transformation. These two preprocessing steps 
result in minimal changes to the ILC map (the maps re- 
main highly correlated with cross-correlations c > 0.95). 
The main effect is the removal of significant outliers in 



^ http://lambda.gsfc.nasa.gov/ 

^ http:/ /space. mit.edu/home/tegmark/wmap. html 
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with 1 < / < Icut^ < m < I are shuffled while the aheady 
randomised phases for the smah scales are preserved. 
Fig. [l] shows a realisation of a second order surrogate 
map after inverse Fourier transformation. Note that the 
Gaussian properties of the remapped ILC map, which are 
given by \aim\^ are ea;ac/^ preserved in all surrogate maps. 
Finally, for calculating higher order statistics the maps 
were degraded to Nside = 256 and residual monopole and 
dipole contributions were subtracted. To compare the 




FIG. 1: ILC map after remapping of the temperatures and 
phases (above). First order (middle) and respective second 
order surrogate (below) for Icut = 20. Note the resemblance 
of the first order surrogate with the ILC map at large scales 



the temperature distribution. To test for scale-dependent 
non-Gaussianities in a model-independent way we pro- 
pose the following two-step procedure. Without loss of 
generality we restrict the description of the method and 
all subsequent analyses to the case of non-Gaussianities 
on large scales. Consider a CMB map T(^,^), where 
T{0^(j)) is Gaussian distributed and calculate its Fourier 
transform. The complex valued Fourier coefflcients a/m, 
o^im = / dO^nT {n)Y^^{n) can be written as aim = 
l^^zmle'"^'^ with (l)im = ^Tct^n{Im{aim)/Re{aim))- The 
linear or Gaussian properties of the underlying random 
field are contained in the absolute values |a^^|, whereas 
all HOCs - if present - are encoded in the phases ^^^n 
and the correlations among them. First, we generate a 
first order surrogate map, in which any phase correlations 
for the scales, which are not of interest (here: the small 
scales), are randomised. This is achieved by a random 
shuffle of the phases for / > /cnt^O < m < where 
hut = 10,15,20,25,30 in this Letter and by performing 
an inverse Fourier transformation (Fig. [T]). Second, N 
{N = 500 for Icut = 20, N = 100 otherwise) realisations 
of second order surrogate maps are generated for the fflst 
order surrogate map, in which the remaining phases (j)im 




FIG. 2: Deviation S as derived from rotated upper hemi- 
spheres for (Jt (above) and (a(rio)) (below) for the WMAP5 
map and Icut = 20. The z-axis of the respective rotated ref- 
erence frame pierces the centre of the respective colour-coded 
pixel. 768 rotated hemispheres, which correspond to num- 
ber of coloured pixels, were considered. (For a more detailed 
description of this visualisation technique see e.g. [14t ll8)). 

two classes of surrogates, we calculate local statistics in 
the spatial domain, namely scaling indices (SIM) as de- 
scribed in Rath et al. [18 . In brief, scaling indices esti- 
mate local scaling properties of a point set P. The spher- 
ical CMB data can be represented as a three-dimensional 
point distribution P = pi = {xi,yi,Zi),i = 1, . . . , Npi^eis 
by transforming the temperature fluctuations into a ra- 
dial jitter. For each point pi the local weighted cu- 
mulative point distribution p is calculated p{pi^r) = 
^N^^.eis ^ ||^_ _^-. ||^ rj.^^ weighted scaling 

indices a{pl^r) are then obtained by calculating the loga- 
rithmic derivative of p{pl,r) with respect to r, a{pl^ r) = 
^^^diogr^^ . For each pixel we calculated scaling indices 
for ten different scales, ri = 0.025,. . . ,rio = 0.25 in 
the notation of [18 . For each scale we calculate the 
mean {{a)) and standard deviation (ctq,) of the scaling 
indices a{pi^r) derived from a set of pixels belonging 
to rotated hemispheres or the full sky. To investigate 
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FIG. 3: Probability density P {a{r lo)) for the surrogates of 
the WMAP5 (blue), T0Hw3 (yellow) and T0Hc3 (red) map 
for the rotated upper and lower hemisphere and I cut = 20. 
The black lines denote the respective first order surrogate. 
The reference frame is chosen such that the difference AS = 
Sup — Slow between the upper and lower hemisphere becomes 
maximal for (a(rio)) regarding the WMAP5 surrogates. 



TABLE I: S/SL Upper Hemisphere 





WMAP5 


T0Hc3 


T0Hw3 




(S/SL) 


(S/SL) 


(S/SL) 




-2.8/ 99.8 


-3.0/>99.8 


-2.9/ 99.8 


(a(rio)) 


3.5 / >99.8 


3.5 / >99.8 


3.6 / >99.8 




5.7 / 99.8 


5.2 / 99.6 


7.0 / >99.8 




3.1 / 99.2 


-0.7 / 74.4 


2.1 / 95.4 




6.1 / > 99.8 


3.6 / 99.0 


6.4 / > 99.8 



Lower Hemisphere 
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FIG. 4: Deviations |6'(r)| for the rotated upper and lower 
hemisphere for (a) (black), cr^ (blue) and a x^-combination 
of (a) and (7^ (red) {Icut = 20, N = 500). The solid (dashed, 
dashed-dotted) lines denote the WMAP5 (T0Hw3, T0Hc3) 
map. The shaded region indicates the 3a significance interval. 
The insets show the results for (a(rio)), (^(rg)) and (^(rs)) 
(solid, dashed, dashed-dotted) as a function of Icut for the 
WMAP5 map (here: N = 100). 



the correlations between the scaling indices and tem- 
perature fluctuations, we also considered the standard 
deviation (ctt) for the mere temperature distribution of 
the respective sky regions. The differences of the two 
classes of surrogates are quantified by the a-normalised 
deviation S{Y) = {Ysurroi - {ysurro2)) / crYsurr-o2i y = 
(jT^{(y)^(Jo,^X^ (surrol: first order surrogate, surro2: sec- 
ond order surrogate) and the significance levels SL = 
1 — p, where p is the fraction of second order surro- 
gates, which have a higher (lower) Y than the first or- 



TABLE II: S/SL Lower Hemisphere 





WMAP5 


T0Hc3 


T0Hw3 




(S/SL) 


(S/SL) 


(S/SL) 


(JT 

(a(rio)> 

Xaa 


2.7/99.8 
-3.9 / >99.8 
7.9 />99.8 
-0.7 / 76.4 
5.8 / 99.8 


2.9/>99.8 
-3.9 / >99.8 
5.4/99.8 
4.4 / 99.6 
6.3 />99.8 


2.8/99.8 
-3.7/>99.8 
7.3/>99.8 

-0.6/67.0 
5.2/>99.8 



der surrogate, denotes diagonal x^-statistics, which 
we obtain by combining {a),aa for a given scale r^, i.e. 

2 

, with Xi — (a),X2 = 



x'(n) = El 



Xjin)-{Xjin)) 

derived from the N realisations of sec- 



aa and {Xj),axj 
ond order surrogates. As scale-independent measure we 
also consider as obtained by summing over the scales 



(Nr = 10), 



one single measure (ji = 1,^2 = 
and the two measures (ji = 1,^2 



X,(r,)-(X,(r,)) 



i) 



for 



(^Xj (r 

1; ji = 2,j2 = 2) 
= 2). Fig. [2] shows 
S(aT) and S ({a(r lo))) derived from pixels belonging to 
the respective upper hemispheres for 768 rotated refer- 
ence frames. Statistically significant signatures for non- 
Gaussianity and ecliptic hemispherical asymmetries be- 
come immediately obvious, whereby these signatures can 
solely be induced by large scale HOCs. Although S((Jt) 
and S({a(rio))) are spatially highly (anti-) correlated 
(c = —0.95), the two effects are nevertheless comple- 
mentary to each other in the sense that a systematically 
lower /higher (Jt would lead to a lower/higher (a(rio)) 
and not to the observed higher/lower value for the first 
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order surrogate map. These systematically shifted scal- 
ing indices are a generic feature present in all three maps 
(Fig. [3|. Although the probability densities P(a(rio)) 
are different due to the smoothing or Wiener-filtering for 
the three maps, the shifts of the first order surrogate 
relative to its second order surrogates can be found in 
all three cases. We also cross-correlated the deviation 
maps shown in Fig. |2] derived from the three input maps 
and always obtained c > 0.98 for the correlation coef- 
ficient. These systematic deviations lead to significant 
detections of non-Gaussianities which are shown in Fig. 
[4] and summarised for Icut — 20 in Tables [mil The most 
significant and most stable results are found for {a) at 
larger radii, where for all three maps none of the 500 sec- 
ond order surrogates had a higher (upper hemisphere) or 
lower (lower hemisphere) value than the respective first 
order surrogate, leading to a significance level SL > 99.8 
% for (a(rio)). Also the combined measure x\a) yields 
deviations S ranging from 5.2 up to 7.9, which represent 
one of the most significant detection of non-Gaussianity 
in the WMAP data to date. We estimated how varying 
Icut values affect the results and found that both the non- 
Gaussianities and asymmetries are detected for all con- 
sidered Icuti where the highest deviations are obtained for 
hut = 20. Although S becomes considerably smaller for 
hut = IO7 we can still detect the non-Gaussianities with 
SL > 99.0 %, which is larger than the results reported 
in [28] {SL = 95 %), where also hut = 10 was used. 
We perfomed the same analyses for the coadded WMAP 
foreground template maps and for simulations using the 
best fit ACDM power spectrum and WMAP-like noise 
and beam properties. We found in none of these cases 
significant signatures as reported above. Details about 
these studies are deferred to a longer forthcoming publi- 



cation. 

In conclusion, we demonstrated the feasibility to gen- 
erate new classes of surrogate data sets preserving the 
power spectrum and partly the information contained in 
the Fourier phases, while all other HOCs are randomised. 
We found significant evidence for both asymmetries and 
non-Gaussianities on large scales in the WMAP data of 
the CMB using scaling indices as test statistics. The 
novel statistical test involving new classes of surrogates 
allows for an unambigous relation of the signatures iden- 
tified in real space with scale-dependent HOCs, which 
are encoded in the respective Fourier phase correlations. 
Our results, which are consistent with previous findings 
[Il[l3[ll[l5l[Il[l8l[28] but also extend to smaller scales 
than those reported in [17] {hut = 3) , [28] {hut = 10) 
and [20J {hut ^ 3), point towards a violation of statis- 
tical isotropy and Gaussianity. Such features would dis- 
favour canonical single-field slow-roll inflation - unless 
there is some undiscovered systematic error in the col- 
lection or reduction of the CMB data or yet unknown 
foreground contributions. Thus, at this stage it is too 
early to claim the detected HOCs as cosmological and 
further tests are required to elucidate the true origin 
of the detected anomalies. Their existence in the three 
maps might, however, be suggestive. 
In either case the proposed statistical method offers an ef- 
ficient tool to develop model-independent tests for scale- 
dependent non-Gaussianities. Due to the generality of 
this technique it can be applied to any signal, for which 
the analysis of scale-dependent HOCs is of interest. 
Many of the results in this paper have been obtained us- 
ing HEALPix [29|. We acknowledge the use of LAMBDA. 
Support for LAMBDA is provided by the NASA Office 
of Space Science. 
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